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The forthcoming LHC measurement of hadron multiplicity in proton-nucleus collisions is a cru- 
cial test of the kt factorization and gluon saturation based models. Here, we provide quantita- 
tive predictions for the pseudorapidity distribution of charged particles produced in minimum bias 
proton-nucleus collisions at the LHC based on the idea of gluon saturation in the color-glass con- 
densate framework. Our formulation gives good descriptions of the LHC and RHIC data for the 
charged-hadron multiplicities in both proton-proton and nucleus- nucleus collisions, and also the 
deep inelastic scattering at HERA at small Bjorken-x. 



I. INTRODUCTION 



The Color Glass Condensate (CGC) formalism is a self-consistent effective QCD theory at high energy (or small 
x) in which one systematically re-sums quantum corrections which are enhanced by large logarithms of \/x and 
also incorporates non-linear high gluon density effects which are important where the physics of gluon saturation is 
dominant, for a review see Ref. [l| and references therein. 

In the CGC approach, the hadron production in proton-nucleus collisions goes in two stages: production of gluons 
' and subsequently the decay of gluon-jet or mini-jet into hadrons, namely hadronization (or fragmentation in the 
standard pQCD language) The first stage of this process is under theoretical control and can be calculated via 
O i' the kt factorization. The jet decay, unfortunately, can be treated mostly phenomenologically. However, one may 
^ , hope that the phenomcnological uncertainties would be reduced to few constants whose values will be extracted from 
the experiment in other reactions and energy. Notice that the inclusive giuon-production in the kt factorization has 
] been proven [2| at the leading log(l/x) approximation for scatterings of a dilute partonic system on a dense one 
CNJ ■ such as p-|-A collisions at high energy-'^. In nucleus-nucleus collisions at around midrapidity, one has to deal with a 
' dense-dense scatterings, therefore the kt factorization and as a consequence all saturation based predictions are less 
reliable. Therefore, the upcoming data from p-l-A run at the LHC is very crucial in order to test the kt factorization 
"J ' and discriminate between various saturation models. 

, Already the first LHC data on hadron multiplicity in nucleus-nucleus collisions call for a theoretical understanding 
• ■ of these data based on QCD. One of the unexpected new feature of the LHC data has been the very different power- 
' law energy behavior of charged hadron multiplicities in A-t-A compared to p-l-p collisions fsl, l6| . Although this is still 
^ ■ open problem, there have been some different approaches to accommodate this feature of data within the saturation 
\ picture 0, Q ■ We will show below that the first-day experimental p+ A data on charged hadron multiplicity can also 
• • . discriminate between these approaches. 

In this letter, we provide quantitative predictions for the pseudorapidity distribution of charged particles produced 
in minimum bias proton-nucleus collisions at the LHC based on the kt factorization within the CGC framework. This 
• approach has been also successful to describe the charged hadron multiplicity in both p-l-p and A-|-A collisions at the 
5^ ^ LHC 0, ll], see also We refer the reader to Refs. d, d, [i3 for the technical details. Here we shortly introduce 

the main formalism in Sec. II and then show our results in Sec HI. 



II. MAIN FORMULATION 



The gluon jet production in A-l-B collisions can be described by /cT-factorization given by 0] , 

da 2a^ 1 



rf^fcr'/'A (a^iifcr) 0s (a;2;PT - fcr) , (1) 



where Cp = {N^ — \)/2Nc with Nc being the number of colors, xi^2 = {pt / \/s)e^'^ , px and y are the transverse- 
momentum and rapidity of the produced gluon jet. (l)'^{xi\kT) denotes the unintegrated gluon density and is the 



^ For the recent theoretical development on this line including the effect of the running coupling corrections to the kt factorization, see 
Ref. g. 
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probability to find a gluon that carries Xi fraction of energy with kx transverse momentum in the projectile A (or 
target B). The unintegrated gluon density is related to the color dipole forward scattering amplitude, 

^2 (x,; kr) = J dH£rTe^^^--^^-VlN'i [xv^tt; b) , (2) 

with notation 

(x,; tt; b) = 2Na {x.^tt; h) - Nl {x,]rT] b) , (3) 

where rj- denotes the dipole transverse size and b is the impact parameter of the scattering. For the value of strong- 
coupling Us we employ the running coupling prescription used in Ref. [l2j . 

The most important ingredient of the single inclusive hadron production cross section which captures the saturation 
dynamics is the fundamental (or adjoint) dipole cross section, the imaginary part of the forward quark anti-quark 
scattering amplitude on a proton or nucleus target Na.v {xj'.rT'-, b). In principle, the dipole cross-section can be 
computed via the JIMWLK/BK evolution equations [ij, 13 • For the recent progress in obtaining the impact- 
parameter dependent solution of the BK equation, see Ref. [l^. It has been shown that the impact-parameter 
dependence of the dipole amplitude is important for phenomenological studies [1, 0, H, [13, , including for describing 
the diffractive data at HERA [l^. Following Refs. [SB Ell ^^e the b-CGC saturation model [l^ which explicitly 
depends on impact-parameter. This model approximately incorporates all known properties of the exact solution to 
the BK equation including the impact-parameter dependence of the scattering amplitude ^] and it also describes 
the small-x data at HERA [16]. For comparison, we will confront the predictions coming from the b-CGC with the 
MCrcBK [3|- The MCrcBK model effectively incorporates the impact-parameter dependence into the BK equation 
solution via the Monte Carlo implementation of the kt factorization by taking into account fluctuations of the position 
of nucleon inside a nucleus [2]. Notice that both the b-CGC and the MCrcBK model have been quite successful to 
describe various experimental data in different reactions. However, we will show here that these two models give 
rather different predictions for the charged hadron multiplicity in p-f A collisions at the LHC. 

The kt factorization has infrared divergence. By introducing a new parameter rrijet as mini-jet mass which mimics 
the pre-hadronization effect, one can also regularize the cross-section. Unfortunately, we do not know how mini-jet 
mass changes with medium and kinematics. Therefore, one may model or chose the value of mini-jet at lower energy 
and hope that it should be still valid at higher energy for different rapidities and centralities. This is the approach 
that all previous saturation based studies have adopted, although, in different models, different values for rrijet has 
been taken [1, [3, [1, [IS [H, [13 ■ Therefore, the uncertainties associated with our freedom in choosing the mini-jet mass 
value should be considered and will be quantified here. In order to take account of the difference between rapidity 
y and the measured pseudo-rapidity 77, we employ the Jacobian transformation between w and 77 [3|. The final-state 
hadronization effects are incorporated via the Local Parton-Hadron Duality principle [s, [3, [1, Iiol . [T^ , namely the 
hadronization is a soft process and cannot change the direction of the emitted radiation^. This introduces only one 
extra parameter which can be absorbed into the over-all normalization of the cross-section. Therefore, we have only 
two unknown parameters in our model, the overall normalization factor and the mini-jet mass which are fixed at lower 
energy at midrapidity for central collisions. Then our results at higher energies, different rapidities and centralities 
can be considered as free-parameter predictions. 



III. DISCUSSION AND PREDICTIONS 



First, in Fig. [TJ we show our description of the existing experimental data in p-l-p inelastic non-singlet diffractive 
(NSD) collisions and predictions for higher energies. Note that our curve at 7 TeV was prediction and it is seen that 



it is in good agreement with the recent LHC data 18|. In Fig. [T] (right), we also show our predictions for charged 



hadron multiplicity distribution in p-l-p collisions at ^/s = 4.4, 14 TeV (s is the center-of-mass energy squared per 
nucleon) . 

In the case of scatterings on a nucleus target, one should also know the atomic number (A) and impact-parameter 
dependence of the saturation scale. In our approach, the saturation scale of proton and nucleus both depend on the 
impact-parameter, see Refs. [1, [l3- However, we assume that the atomic number or A (size of nuclei) dependence of 



^ One should note that the main contribution of the kt factorization for the multiplicity comes from small pj- < 1.5 GeV where the 
fragmentation functions are not reliable. We recall that it is still an open problem how to incorporate the fragmentation processes into 
the CGC/saturation formalism. 
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FIG. 1: Pseudo-rapidity distribution of charged particles production in p+p (right panel) and A+A (left panel) collisions at 
various energies. Right: the theoretical curves from top to bottom are at ^/s = 14, 7 and 4.4 TeV. The band in the right panel 
indicates about 2% theoretical error. The total theoretical uncertainties is less than 6% mainly due to our freedom in fixing 
the overall normalization and mini-jet mass with the experimental data at lower energy, see Ref. 01 . The experimental data 
are from [1,1, [Til. 



the saturation scale is approximately factorizable from energy or x for a large nuclei [8|, [12| . This assumption may be 
simply justified from the experimental observation in A+A collisions that the centrality dependence of the multiplicity 
at the LHC is very similar to RHIC up to an overall normalization Moreover, the pseudorapidity distribution 
of multiplicity at a fixed energy but different centralities approximately have similar shape upto a normalization 
factor |8| . These experimental facts indicate that the energy and centrality (or A) dependence of multiplicity, and as 
consequence the saturation scale, can be approximately scaled out, in accordance with our assumption. 

As we already pointed out, one of the most interesting new feature of the LHC multiplicity data is the fact 
that the energy growth of multiplicities in A+A collisions is different from p+p ones [1, [a]. Notice that the kt 
factorization at the leading log approximation accounts for most of hadron multiplicity while still some contributions 
due to gluon decay cascade in the so-called MLLA (Modified Leading Logarithmic Approximation) [l^ kinematics 
may be missing [8]. The latter effect requires the inclusion of higher order corrections and has been traditionally 
studied within a different resummation scheme in which one systematically incorporates next-to-leading logarithmic 
corrections, namely single and double- logarithmic effects in the development of parton cascades [l^. In Ref. [sj we 
extracted the energy-dependence of the gluon-jet decay cascade from e+e" annihilation data and we showed that the 
energy-dependence of about s" °'^^ due to the enhanced gluon-decay cascade is exactly what explains the different 
power-law energy-dependence of hadron multiplicities in A+A compared to p+p collisions at the LHC. This effect is 
more important for A+A collisions at high energy where the saturation scale is larger and consequently the average 
transverse momentum of the jet becomes about or bigger than 1 GeV. The MLLA gluon decay cascade gives rise to 
an extra contribution about 20 — 25% to the multiplicity in A+A collisions at the LHC. This picture is fully consistent 
with all existing experimental data including RHIC data We will show below that this picture can be further 
tested in p+A run at the LHC. 

In Fig. [1] (left), we show the charged hadron pseudorapidity distribution at the LHC in A+A central (^ — 5%) 
collisions at various energies. The effect of the gluon-decay cascade are incorporated in our results. Again, we have 
only two parameters which was fixed with the lower energy data points at the RHIC. Although, the midrapidity 
experimental data point at i/s = 2.76 TeV for A+A — 5% central collisions was available [3| at the time that we 
published our results [1] , the pseudorapidity distribution was prediction and its agreement with the recently released 
data from the CMS collaboration [g^ is striking. 

Having described both p+p and A+A multiplicities experimental data for large range of energies/rapidities, the 
multiplicity in p+A collisions can be considered as an important check of consistence. The normalization again 
(similar to Fig. [Ij is fixed with RHIC data. Unfortunately the error bars of RHIC data is rather large (see Fig. [5]) and 
this induces uncertainties upto ~ 15% in our results. It was shown in Refs. ^, i8(] that assuming a different mini-jet 
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FIG. 2: Pseudo-rapidity distribution of the charged particles production in p+A minimum (Mini) bias colhsions at the 
LHC \/s = 4.4 TeV. The theoretical curves labeled by CGC-Rl and MCrcBK are based on the leading log fct-factorization 
formalism but two different saturation models. The curve labeled by CGC-R2 is based on the modified fcj -factorization formalism 
incorporating the gluon-decay cascade effects from the Modified Leading Logarithmic Approximation scheme, see the text for 
the details. The experimental data are from (20(). 

rrijet for the cases of A+A and p+p collisions will be in favor of data. Here, in order to reduce the model dependence 
associated with the choice of mini-jet mass, we take a fixed value for mj^.^ for all rapidities, in the both fragmentation 
region of nucleus and proton. We estimated that uncertainties due to various value for rrijet brings about 10% errors 
in our calculation. We focus only on a pseudorapidity region where the effect of valence quarks can be ignored. In 
Fig. [21 we show our prediction for the pseudorapidity distribution of charged hadron multiplicity at the LHC energy 
^/s = 4.4 TeV in minimum-bias p-|-Pb collisions. We show our results both based on the kt factorization (labeled by 
CGC-Rl) and also the modified kt factorization incorporating the effect of the MLLA gluon-decay cascade (labeled by 
CGC-R2). Notice that the contribution of the MLLA gluon-decay cascade may be also important in p+A collisions 
at the LHC energy since the average saturation scale of the system can be bigger than 1 GeV. This effect enhances 
the multiplicity about 20%, in accordance with our results for A+A collisions Q. 

Notice that in the case of p+p and A+A collisions where the interacting system is symmetric in the rapidity, the 
corresponding multiplicity distribution is also symmetric (see Fig. [1}. The position of two peaks in pseudorapidity 
corresponds to the fragmentation region of the projectile and the target. It has been shown that the peak in the 
multiplicity distribution at forward and backward rapidity is enhanced at higher energies due to the fact that the 
saturation scale increases with energy (and density) 0, H, O, [l^, [13, HH . The exact shape of multiplicity distribution 
(position of the peaks, the local minimum and the width of the peak) depends on the unintegrated gluon density 
and the Jacobian transformation which relates the rapidity and pseudorapidity. In the case of p+A collisions, the 
interacting system is not symmetric namely the saturation scale in the nuclear fragmentation region is bigger than the 
projectile (proton) side and consequently the multiplicity distribution in pseudorapidity is not symmetric. In Fig. [51 
we see rather strong forward-backward asymmetry for p+Pb compared to the symmetric collisions. The shape and 
position of maximum in backward and forward rapidity distribution are intricately related to the saturation scale in 
nucleus and proton. A bigger saturation scale in the nucleus washes away the maximum of the multiplicity distribution 
in the proton fragmentation region. 

For comparison, in Fig. [2l we also show the predictions coming from the MCrcBK of Ref. [7] . Assuming that both 
predictions have about 15% errors due to the uncertainties associated with fixing the normalization with the RHIC 
data, it is seen that two approaches give rather different predictions for the multiplicity distribution in p+A collisions. 
Using the kt factorization at the leading log approximation, our prediction for charged hadron multiplicity at ^/s = 4.4 
TeV at midrapidity = in p+A mini-bias collisions is dN/dr] = 19.67 ± 1.5 (CGC-Rl) while prediction coming from 
the MCrcBK approach is dN/drj = 16.71 ± 1.3 (MCrcBK). The main difference between these two approaches is 
due to different employed saturation model. On the other hand, including the MLLA gluon decay cascade effect into 
the kt factorization increases our multiplicity results and at midrapidity we have dN/d-q = 21.88± 1.7 (CGC-R2), see 
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Fig.H 

To conclude, in this letter, we provided quantitative predictions for charged hadron multiplicity in p+Pb collisions 
at the LHC. We showed that the LHC p+Pb data on the charged particles pseudo-rapidity multiplicity distribution 
can discriminate between saturation models and also examine the kt factorization at an unprecedented level. In 
general, the measurement of inclusive hadron production in p+A collisions at the LHC can be a very good probe of 
the small-x physics, see also Ref. (22l |. 
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